PACS numbers: 36.40. Qv, 71.15.Nc, 71.15.Pd. 

Metastable spin-polarized carbon clusters and their ensembles 

A. A. Izmalkov, L. A. Openov 

Moscow Engineering Physics Institute (State University), 
Kashirskoe sh. 31, Moscow 115409, Russia 

Results of tight-binding calculations on metastable spin-polarized 
carbon clusters Cs (prismanes) are presented. It is shown that those 
clusters can form ensembles due to intercluster bonding. The binding 
energy of a given metastable configuration decreases monotonously with 
the total spin S, while the activation energy for the decay of a metastable 
state has distinct maxima as a function of S. For specific values of S, 
the dynamical stability of spin-polarized ensembles of prismanes appears 
to be higher than that of an isolated prismane, pointing to a possibility 
of existence of "cluster matter" composed of spin-polarized clusters Cs. 

1. Introduction 

Carbon is a chemical element with a great number of unique properties. In par- 
ticular, carbon is known by its polymorphism. Due to its ability to exist in different 
valence states, carbon can form a rich variety of compounds, including diamond, 
graphite, carbyne, etc. In the field of nanostructures, cage-like fuUerenes (Ceo etc. 

attract much attention due to their both practical and fundamental importance. 
Recently, the cage-like cluster C20 was synthesized 0. This cluster is the smallest 
possible fullerenes cage consisting of 12 pentagons. First-principle calculations 
point to a possibility of synthesis of condensed C20 phases with three-dimensional 
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covalent networks. 

In recent numerical simulations of small carbon clusters, a possible existence of 
another cage-like cluster Cg (prismane) was suggested theoretically [Q. Prismane 
was predicted to have the shape of a six-atom triangular prism with two excess 
atoms above and below its bases, see Fig.l. There is a good reason to believe that 
prismane Cs is the smallest possible three-dimensional carbon cluster. Prismane is 
the metastable cluster. Its binding energy Eb (defined as Ef,{N) = NE{1) — E{N), 
where E{N) is the total energy of an A^-atom cluster, E{1) is the energy of an 
isolated carbon atom) was calculated to be 5.1 eV/atom 0], i.e., 0.45 eV/atom 
lower than the binding energy of the stable one-dimensional eight-atom cluster and 
2.3 eV/atom lower than the binding energy of the bulk graphite or diamond. It 
is essential, however, that the calculated value of the lowest energy barrier that 
prevents the prismane from decay into the stable or another metastable atomic 
configurations with greater values of Efy appears to be rather high, about 0.44 eV 
0, resulting in a rather long lifetime of the prismane even at room temperature. 
Hence, it is likely that this cluster may be observed experimentally. 

Since the prismane had been predicted the bonding between prismanes was 
studied and several configurations in which prismanes form metastable ensembles 
were found [0, |3, |]. In those ensembles, the inter-prismane bonds are rather strong, 
so that one could expect prismanes Cg to form a covalent solid. However, the 
stability (and hence, the lifetime) of the prismane-based ensembles was shown to 
decrease strongly with the number of prismanes in the ensemble. 

Until now, theoretical studies of prismanes and their ensembles were restricted to 
the non-polarized case, i.e., the total spin 5* of the systems under consideration was 
assumed to be zero 0, ||, ^ 0, H. Meanwhile, the polarization of small species of 
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substance (including molecules and clusters) can increase the stability of the system 
and/or result in the formation of new metastable configurations. For example, two 
helium atoms, being in the ground state with S = each, do not form the covalently 
bound molecule He2, while the existence of metastable long-lived molecule Heg in 
the excited state with S = 1 has been demonstrated both theoretically ^ and 
experimentally |10[- In general, one can expect that the binding energy of the system 



composed of nonmagnetic atoms decreases with S, and hence spin-polarized clusters 
and their ensembles can be considered as candidates to the so called high energy 
density materials (HEDM). The question is whether the lifetimes of such systems 
are long enough for they could be synthesized, investigated and used in practice. So, 
the study of metastable spin-polarized atomic configurations is of great fundamental 
interest and practical importance. 

In this paper, we extend our previous studies to the case of spin-polarized pris- 
manes Cg and their ensembles. We show that metastable prismane-like carbon 
clusters with S can exist and form ensembles due to intercluster bonding. The 
binding energies and the activation energies for the decay of metastable states are 
calculated as functions of S for an isolated prismane and different arrangements of 
prismanes in the ensembles. 

2. Computational details 

The search for metastable configurations of prismane-like carbon clusters and their 
ensembles with different spins S was carried out by means of structural relaxation 
technique. Starting with a particular atomic arrangement, all atoms were allowed 
to displace step by step in the directions of interatomic forces, atomic velocities 
being terminated at each step. The forces were determined via calculation of the 
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total energy of the system (or, equivalently, the binding energy) making use of 
a transferable tight-binding potential [IT, T2\ that had been proven to reproduce 



accurately the energy-versus-volume diagram of carbon polytypes and to give a 



good description of both small clusters and bulk structures of carbon [pT| , |l^. This 
potential has been used for simulations of metastable prismane-like structures with 
5* = in our previous works ^ ||, |[. The numerical technique employed is 
not as sophisticated and time-consuming as ab initio ones and allows for molecular 
dynamics simulations of rather large systems for a relatively long time. 

Depending on the choice of initial configuration, the resulting atomic arrangement 
corresponds to either global or local minimum of the total energy as a function of 
atomic coordinates for a given number of atoms in the system (stable and metastable 
states, respectively). We were interested in metastable states separated from the 
stable and/or other low-lying metastable configurations by the energy barriers. An 
important characteristic of a metastable configuration is the activation energy Ea 
for its decay. For a given temperature T, the value of Ea determines the probability 
of the system decay in a unit of time, 

W oc exp {-EJUbT) , (1) 

and hence the lifetime r of a particular metastable state is 

T (XTq exp {Ea/ksT) , (2) 

where fc^ is the Boltzmann constant, tq is a characteristic time of the order of 
an oscillation period of the system (~ 10^^^ s for the Cs clusters and their small 
ensembles) . 

The lifetimes r at different temperatures were calculated directly by molecular 
dynamics simulation of thermally isolated systems. Since the cluster decay is prob- 



abilistic in nature, we accumulated a great statistics (30 100 trials) in order to get 
a reliable estimate of the value of Ea through fitting the results obtained to Eq. (|^). 

3. Results and discussion 

To find three-dimensional spin-polarized prismane-like configurations of eight car- 
bon atoms, we made use of the non-polarized prismane structure as an initial 
atomic arrangement in the structural relaxation procedure. We have found three- 
dimensional metastable Cg configurations for 5 = 1 7, while for 5* > 7 all relaxed 
structures appeared to be either one- or two-dimensional. Metastable configurations 
can be subdivided into "nondistorted" {S = 0, 4, 7) and "distorted" {S = 1, 2, 3, 5, 6) 
ones. The structure of nondistorted prismanes is characterized by three different 
bond lengths a, b, c (symmetrical configuration), see Fig.l. In distorted prismanes, 
the bond lengths are slightly scattered around those three characteristic lengths, see 
Fig. 2. From Fig.2 one can see that the values of a, b, and c (or their average values 
for distorted configurations) increase with S, so that the cluster volume increases 
with S. 

The binding energy Eh of spin-polarized prismanes strongly decreases with S 
from Efy = 5.1 eV/atom at S" = down to E^ = 2.1 eV/atom at S* = 7, see the 
lower curve in Fig.3. Hence, the energy Eacc accumulated in metastable prismanes 
increases considerably and exceeds 5 eV/atom at 5* = 7 (here Eacc = E^^"'^^^^'^ — Eb, 
where E^raphtte _ eV/atom is the binding energy of graphite or diamond, the 
most stable carbon structures that are characterized by the maximum value of Eb 
per atom). Such an energy could be released upon fusion of a large number of 
prismanes into the bulk specimen with the graphite structure, accompanying by the 
spin-depolarization transition (i.e., the reduction of the total spin down to S* = 0). 
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Fig. 4 shows the dependence of activation energy Ea for the decay of metastable 
spin-polarized prismane Cg versus the total spin S. One can see that the value of 
Ea is maximum, Ea = (0.84 ± 0.05) eV, for the non-polarized prismane Cg with 
S = 0. An estimate of the lifetime r making use of Eq. (^) gives r ~ 1 s at T = 273 
K (see also Ref. [^), the value of r being increased exponentially with decreasing 
T. For spin-polarized Cg configurations, the value of ii^a is 4 12 times lower than 
that for S* = 0. Note that the dependence of Ea on 5" in the range S = 1 ^ 7 is 
nonmonotonous and has local maxima at 5 = 4 and S = 7. 

In attempts to get a physical insight into the nonmonotonous dependence of Ea 
versus S, we have calculated the difference A = Em mo ~ Ehomo between the 
energies of the lowest unoccupied and the highest occupied molecular orbitals as a 
function of 5*. The results are shown in Fig.5. One can see that the curve A(5') has 
maxima at S = 0, 4, and 7, in accordance with higher values of Ea for those values 
of S, see Fig.4. 

Now we turn to small ensembles composed of spin-polarized prismanes Cg. We 
have studied various arrangements of two prismanes Cg with respect to each other 
and found four different metastable configurations (Cg)2. They are shown in Fig.6. 
In configuration A, the intercluster bonding occurs via "top atoms", so that there 
is a single inter-cluster bond. In configurations B and C, there are two inter-cluster 
bonds between "edge atoms", while in configuration D, there are four inter-cluster 
bonds (the rectangular facets of the prismanes are parallel). Note that prismanes in 
the ensembles (Cg)2 preserve their overall original shape, they do not merge into a 
new cluster Cig and survive as the "building blocks" of the ensembles. 

For either of four configurations (Cg)2 shown in Fig.6, we have restricted our 
studies to the range 5 = 10. We have found that the metastable configuration 
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A exists at all values oi S = 10, while the configurations B, C, and D are 
metastable, respectively, at = and 7; 6 and 7; 8 and 9 only. The binding energy 
Ef, of (03)2 ensembles monotonously decreases with S from Ei, = 5.2 eV/atom at 
S — down to Eb — 3.4 eV/atom at S — 10, see the middle curve in Fig. 3. The 
value of Ef, at a given S is almost the same in different metastable configurations 
(03)2, if any 

Although at 5" 7^ the curve Eb{S) for the ensembles (C8)2 lies well above that 
for the prismane Cg, see Fig. 3, a close inspection of those curves shows that the 
binding energy Ei,{2S) of an ensemble (03)2 composed of two prismanes Cs with the 
spin S each is almost the same as the binding energy Ei,{S) of an isolated prismane 
Cg with the spin S, being higher by just (0.1 0.2) eV/atom for all values of S 
considered. Hence, the energy Eacc — E^^"'^^^*''^ — Eh accumulated in two metastable 
prismanes Cg changes insignificantly upon formation of the ensemble (Cg)2. A small 
increase in Ei, by (0.1 0.2) eV/atom is due to appearance of new covalent bonds 
between two prismanes Cg and, as a consequence, the decrease of the total energy 
of the system. The bond energy is about 1 eV/bond. 

In contrast to the prismane Cg, the activation energy £'„ for the decay of the 
ensemble (Cg)2 in configuration A as a function of S has a broad maximum, Ea — 
(0.35 ± 0.05) eV, at 5" = 6 ^ 8, while the value of Ea for the non-polarized ensemble 
(C8)2 with 5" = is very low, Ea ~ 0.03 eV, see Fig.7. (We recall that in the 
prismane Cg, the maximum of the Ea{S) curve is at = 0). For other metastable 
(C8)2 configurations shown in Fig. 6, the value of Ea equals to (0.83 ± 0.12) eV and 
(0.12 ±0.03) eV for configuration B with ,5 = and 7 respectively; (0.03 ±0.005) eV 
and (0.29 ± 0.03) eV for configuration C with S = 6 and 7 respectively; (0.1 ± 0.03) 
eV and (0.07 ± 0.01) eV for configuration D with S — S and 9 respectively. It is 
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worth noting that several spin-polarized configurations (03)2 with specific values of 
S appear to be more stable than constituting spin-polarized prismanes Cg, c.f. Figs. 

4 and 7. 

We have also studied the metastability of quasi-one-dimensional spin-polarized 
ensembles (Cgja formed upon attachment of a third prismane Gs to the ensemble 
(C8)2 in configuration A, see Fig.8. We have restricted ourselves to the range S — 
16. We have found that the metastable configuration (C8)3 exists at all values 
of S in this range. The binding energy of the (Cgjs ensemble monotonously 
decreases with S from — 5.3 eV/atom at 5" = down to Ei, — 3.3 eV/atom at 

5 = 16, see the upper curve in Fig. 3. From Fig. 3 one can find that the binding 
energy Ei,{3S) of an ensemble (03)3 composed of three prismanes Cg with the spin 
S each exceeds the binding energy Eh{S) of an isolated prismane Cg with the spin 
S hy ^ 0.2 eV/atom for all values of S in the range 5" = 1 5. Hence, an increase 
in the number m of prismanes in a metastable spin-polarized ensemble (Cg)^ leaves 
the energy Eacc = Ef^^^^^^ — Ei, accumulated in the ensemble almost unchanged. 

The shape of the curve Ea{S) for the ensemble (Cg)3 is similar to that for the 
ensemble (Cg)2, see Fig.7. Note, however, that the maximum of Ea{S) for the 
ensemble (Cg)3 is shifted to greater values of 5" = 9, 10 and is higher by ^ 0.1 eV 
than in the case of the ensemble (Cg)2. So, the stability (and hence, the hfetime) of 
the system increases with the number of prismanes in the ensemble. 

Fig. 9 shows the activation energy E^ versus the difference A = E^umo ~ Ehomo 
between the energies of the lowest unoccupied and the highest occupied molecular 
orbitals for all metastable configurations considered in this work. One can see the 
overall increase of Ea with A. Such a correlation between Ea and A indicates 
that the nonmonotonous dependence of Ea on S stems from the nonmonotonous 
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dependence of A on S. 
4. Conclusions 

We have numerically examined a possibility of existence of a spin-polarized cage- 
like metastable carbon cluster prismane Cg and small ensembles composed of those 
clusters. We have shown by molecular dynamics simulations that the activation 
energy £'„ for the decay of a metastable configuration is a nonmonotonous function 
of the total spin S of the system and has maxima at certain values of S. The 
most stable structures are characterized by the large energy gap between the lowest 
unoccupied and the highest occupied molecular orbitals. The maximum value of Ea 
increases with the number of prismanes in the ensemble, pointing to a possibility of 
existence of "cluster matter" composed of spin-polarized clusters Cg. The metastable 
ensembles of prismanes Cg accumulate a great amount of energy (several eV/atom) 
that may be released upon fusion of a large number of prismanes into the bulk 
specimen. Hence, such ensembles can be considered as candidates to the high energy 
density materials (HEDM). The role of spin- flip processes for the stability of spin- 
polarized cluster configurations needs further consideration. 
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Figure captions 

Fig.l. Metastable prismane Cs- For 5" = 0, 4, and 7 the prismane structure is char- 
acterized by three different bond lengths a=AE=BF=CG, 6=AD=BD=CD=EH=FH 
and c=AB=AC=BC=EF=EG=FG (symmetrical configuration). For other values 
of S, the prismane is shghtly distorted, see text and Fig.2 for details, a — 1.28 A, 
b = 1.47 A, c = 2.31 A for 5 = 0. 

Fig.2. Bond lengths a (circles), b (triangles), and c (squares) of metastable pris- 
mane configurations versus the total spin S, see Fig.l for the definition of a, b, c. For 
the distorted configurations with S — 1, 2, 3, 5, 6 there are more than three different 
bond lengths, i.e., some of equalities AE=BF=CG, AD=BD=CD=EH=FH=GH, 
AB=AC=BC=EF=EG=FG are violated. 

Fig. 3. Binding energy of the metastable prismane Cg (closed circles), see 
Fig.l, the metastable ensembles {Cs)2 in configurations A (closed triangles), B (open 
circles), C (open squares), D (open triangles), see Fig. 6, and the metastable ensemble 
{Cs)3 (closed squares), see Fig. 8 versus the total spin S of the system. Lines are the 
guides to the eye. 

Fig.4. Activation energy Ea for the decay of the metastable prismane Cg versus 
the total spin S. Line is the guide to the eye. 

Fig. 5. The difference A = Elumo — Ehomo between the energies of the lowest 
unoccupied and the highest occupied molecular orbitals of the metastable prismane 
Cg versus the total spin S. Line is the guide to the eye. 

Fig. 6. Four different configurations (A, B, C, D) of metastable ensembles (C8)2- 

Fig. 7. Activation energy Ea for the decay of the metastable ensemble {C^)2 in 
configuration A (circles), see Fig.6a, and quasi-one-dimensional ensemble {C^)z (tri- 
angles), see Fig.8, versus the total spin S of the system. Lines are the guides to the 
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eye. Two values of Ea of the ensemble (C8)2 at 5 = 1 correspond to two basically 
different paths of the ensemble decay. 

Fig. 8. Quasi-one-dimensional metastable ensemble (08)3. 

Fig. 9. Activation energy Ea for the decay of metastable spin-polarized config- 
urations versus the difference A = Elumo — Ehomo between the energies of the 
lowest unoccupied and the highest occupied molecular orbitals for the metastable 
prismane Cs (closed circles), the metastable ensembles (03)2 in configurations A 
(closed triangles), B (open circles), C (open squares), D (open triangles), and the 
metastable ensemble (03)3 (closed squares). 
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